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Abstract. Self- assembled patterns obtained from a drying nanosphere suspension are investigated by com- 
puter simulations and simple experiments. Motivated by the earlier experimental results of Sasaki & Hane 
and Schope, we confirm that more ordered triangular lattice structures can be obtained whenever a mod- 
erate intensity random shaking is applied on the drying system. Computer simulations are realized on an 
improved version of a recently elaborated Burridge-Knopoff type model. Experiments are made following 
the setup of Sasaki and Hane, using ultrasonic radiation as source for controlled shaking. 



PACS. 81.16.Rf Nanoscale pattern formation - 81.07.-b 
and characterization 



1 Introduction 

In order to use less material, to obtain a more effective 
storage of information or easier access to tight and narrow 
spaces, modern areas of engineering needs structures or 
devices as small as possible. Nanostructures became ideal 
candidates in this sense. Thanks to the efforts of nano- 
chemists, nowadays various nanoparticles nearly monodis- 
perse in terms of size, shape, internal structure, and sur- 
face chemistry, can be produced through reliable and stan- 
dard manufacturing processes. These nanoparticles can be 
nanospheres, nano-tubes or colloids, and can be used as 
building blocks for engineering more complex structures. 
Human assisted construction is however very complicated 
and unproductive. Designing thus technologies where the 
building blocks self- assemble is of primary importance. In 
the present contribution a self-assembling system which 
is widely used in NanoSphere Lithography (NSL) [1,2,3, 
4,5,6,7,8J is studied. Inspired by the earlier experimental 
work of Sasaki & Hane |9j and Schope [10J, here it is con- 
firmed by computer simulations and independent experi- 
ments that an extra shaking mechanism imposed on the 
drying system is benefic for engineering structures with 
better practical properties. 

The studied system is a monodisperse polystyrene na- 
nosphere suspension which is dried on a previously pre- 
pared silica substrate. Self-assembly is governed by sev- 
eral forces: (i) capillarity forces, (ii) electrostatic repulsion 
forces between the slightly negatively charged nanospheres 
(iii) hard-core type repulsion forces between the nanospheres 
and (iv) pinning forces acting on the nanospheres. Many 
previous studies have shown that after the drying pro- 
cess is completed the nanospheres self-organize in com- 
plex patterns. The case when monolayer nanosphere pat- 
terns are obtained is of special practical importance. These 
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Fig. 1. Characteristic self-organized patterns obtained after 
the drying process is ended. The figures from the bottom line 
illustrates the triangular and square crystallization phases. Re- 
sults from our own experiments, for experimental details see 



monolayer structures are useful as deposition masks for 
NanoSphere Lithography (NSL). NSL is nowadays recog- 
nized as a powerful fabrication technique to inexpensively 
produce nanoparticle arrays with controlled shape, size, 
and interparticle spacing [SHU]. 
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Ideally, it is expected that the nanospheres will order 
in the compact triangular lattice structure (Figure [Tfc). 
However, apart of this ideal lattice topology (desirable for 
most of the practical purposes), many other structures 
are formed. Usually the patterns presents many disloca- 
tions, voids or clusters (Fig. [lji,b). Sometimes instead of 
the triangular symmetry, the less compact square topol- 
ogy is selected (Fig. [T]i). A fundamental goal for further 
progress in NSL is the development of experimental proto- 
cols to control the interactions, and thereby the ordering 
of nanoparticles on the solid substrates [6|fT2] . Recently 
[T2] , we proposed a Burridge-Knopoff type model [13] for 
understanding the pattern formation mechanism in this 
system. The model proved to be successful in reproducing 
the observed monolayer patterns, and clarified the influ- 
ence of some experimentally controllable parameters. By 
further improving this model it is straightforward to show, 
that applying an extra random force on the system will 
result in more ordered triangular nanosphere structures. 
One can explain thus the interesting experimental results 
of Sasaki & Hane [9] and Schope [10], and get further ev- 
idences that our Burridge-Knopoff model works fine for 
describing this self-assembling phenomena. 

The present paper will investigate thus by large-scale 
computer simulations this phenomena and present also 
further experimental evidences in support of Sasaki & 
Hane's and Schopfe's results. First our Burridge-Knopoff 
type model [12] is presented and improved. Then simula- 
tions with an extra random or periodic shaking force are 
performed. Finally, as an experimental exercise the ran- 
dom shaking is realized by an ultrasound bath in contact 
with the drying nanosphere system. Both the simulations 
and experiments clearly confirms that the extra shaking 
is useful for engineering more ordered structures. 

2 Computer simulations of the patterns 

Pattern formation in the drying nanosphere suspension 
was recently [12] successfully modeled by a relatively sim- 
ple Burridge-Knopoff type model [T3^ [T4]. The model is 
rather similar with the spring-block stick-slip model used 
for describing fragmentation structures obtained in drying 
granular materials in contact with a frictional substrate 
[T5l fT6]. The new feature of the model presented in [12] is 
that a predefined lattice structure is not imposed anymore. 
The model is two-dimensional; its main elements are disks 
which can move on a frictional substrate and springs con- 
necting them (Figure EJ). Disks, all with the same radius R, 
model the nanospheres. The elastic springs reproduces the 
resultant of all the forces acting between the nanospheres 
for small separation. These forces are the lateral capillarity 
forces and the electrostatic repulsion between the slightly 
negatively charged nanospheres. Springs have similar elas- 
tic constants fc, and their length is defined as the distance 
between the perimeters of the connected disks. There is 
an additional almost hard-core type repulsion Fj which 
forbids disks to interpenetrate. The friction (pinning) be- 
tween disks and surface equilibrates a net force less than 
Ff. Whenever the total force acting on a disk exceeds Ff, 
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Fig. 2. Basic elements of the spring-block stick-slip model. 



the disk slips with an over-damped motion. The tension in 
each spring is proportional with the length of the spring 
Fk = fcZ, and has a breaking threshold F&. 

It worth mentioning here that the approach based on 
the modified Burridge-Knopoff model seems realistic only 
for dense nanosphere systems where the average separa- 
tion between the nanospheres is less than their radius. In 
this limit one can find a plausible argument for modeling 
the resultant forces between the nanospheres with har- 
monic forces. Before proceeding further with the descrip- 
tion of the model we briefly present this argumentation. 

Lateral capillarity forces have been studied in much 
detail both for micrometer size and nanometer size par- 
ticles. Theoretical results [T7lfT8l fT9] predict that for rela- 
tively large separation between two spheres partially im- 
mersed in a liquid one would expect an attractive force 
which is monotonically decaying with separation follow- 
ing a 1/r trend (r the distance between the centers of the 
two spheres). Experimental results confirm this theoret- 
ical prediction for micrometer-sized particles [20,21,22 . 
For nanometer sized particles (our case) experiments were 
made only for the case when the two spheres are confined 
in a liquid film [23] . In this case the measured force had a 
non-monotonic variation with the separation of the parti- 
cles. For small distances it increases almost linearly, than 
reaches a maximum after which decays roughly as 1/r 
with interparticle separation. In our case the polystyrene 
nanonspheres are slightly negatively charged, and there is 
thus an additional Coulombian repulsion between them. 
Assuming that the lateral capillarity force decay as 1/r 
and the Coulombian repulsion decay as 1/r 2 one imme- 
diately gets that the resultant force has a non-monotonic 
variation and for small interparticle separation the resul- 
tant force should increases with r (Figure EJ). The easiest 
way to model this resultant force for small interparticle 
separation is by assuming a linear variation with r. This 
trend is also in agreement with the lateral capillarity forces 
measured between nanospheres confined in a liquid film 
[23] . The above heuristic argument is far from being com- 
plete and lacks the numerical prefactors for the strength of 
the Coulombian repulsion and capillarity attraction forces. 
It serves only as a first motivation for the application of 
the simple spring-block approach. In case of loose systems 
(where the interparticle separation is larger than the di- 
ameter of the nanospheres) one should consider also forces 
decaying as 1/r and the relevance of the spring-block ap- 
proach is questionable. It is interesting however, that even 
in this limit (see the pictures in [12]) the simple approach 
based on the Burridge-Knopoff type model works reason- 



F. Jarai-Szabo et al.: Shake induced order in nanosphere systems 



3 




Fig. 3. Forces acting between the nanospheres. On the image 
from the right-side the resultant force is plotted. Positive force 
represents repulsion and negative forces attraction. 



ably well. In all simulations considered in the present work 
we will restrict our self on dense systems with interparti- 
cle separation much smaller than the particles diameter. 
In this case the harmonic force approach (modeling the re- 
sultant force with springs) is thus heuristically motivated 
and leads to realistic structures. 

In the simulation initially disks are randomly distributed 
and connected by a network of springs with elastic con- 
stants ki n i. (Figure [2}d). We put springs between pairs of 
disks, for which the centers can be connected without in- 
tersecting another sphere (this condition will be referred 
later as the geometric condition). An initially pre-stressed 
spring-block network is thus constructed. During each sim- 
ulation step the spring constant is fixed and the system 
relaxes to an equilibrium configuration where the tension 
in each existing spring is lower than the breaking threshold 
Fb, and the total net force acting on each disk is lower in 
magnitude than the slipping threshold Ff . This relaxation 
is realized through several steps: 



1. For all springs the tension 



is compared with 



If 



r k 



2. 



k ~r j J acting on 



> F^, the spring is broken and taken away 
from the system. 

The total forces F\ = Y, P (d ip Ff + F ip 

disks are calculated (the sum is over all the other disks 
p, dip is 1 if the disks are connected by a spring and 
otherwise, the subscripts k and j denotes elastic forces 
from springs and hard-core type repulsion forces be- 
tween disks, respectively). 

3. Each disk is analyzed. If the magnitude of the total 
force \F\\ acting on a disk is bigger than Ff, then 
the disk will slip with an over-damped motion gov- 
erned by viscosity rj, and its position will be changed 
by dr = F\ dt/rj. The repulsive hard-core potential for- 
bids the spheres to slide on each other and the presence 
of viscosity eliminates unrealistic oscillations. 

4. During the motion of a disk it can happen that another 
spring is intersected. This intersected spring will brake 
and will be taken away from the system. 

5. After all disks have been visited in a random order 
and their possible motions done, the springs that ful- 
fill the considered geometrical condition and for which 



the tension is lower than the breaking threshold are 
redone. By this effect the rearrangement of water be- 
tween nanospheres is modeled. 

This concludes one relaxation step. The relaxation is 
continued until a relaxation step is finished without having 
any spring breaking or disk slipping event. After the relax- 
ation is done, we proceed to the next simulation step and 
increase all spring constants by an amount dk. This step 
models the phenomenon that water level of the continu- 
ous film decreases due to evaporation and the meniscus 
accounting for the capillarity forces gets more accentu- 
ated. The system is relaxed for the new spring-constant 
value, and the spring constant is increased again, until 
all springs are broken or a stable limiting configuration is 
reached. 

Several type of boundary conditions can be imposed. 
One possibility is the free boundary condition which can 
be realized in a simple manner by positioning initially the 
disks inside a circle to minimize the effect of edges. An- 
other possibility is to consider fixed boundary conditions. 
This can be realized for example by positioning again the 
disks inside a circle and considering a chain of fixed disks 
on the chosen circle. These fixed disks are then connected 
with the geometrically allowed springs to other disks. One 
can also consider periodic boundary condition and posi- 
tion initially the disks inside a rectangle. 

The above sequence of events can be implemented on 
computers and relatively big systems with over 10000 of 
disks can be simulated in reasonable computational time. 
The model, as described above, has several parameters: 
the value of the static friction force Ff, the value of the 
breaking threshold F5 of springs, the initial value of spring 
constants ki n% , the spring constant increasing step dk, the 
value of viscosity 77, the parameters of the Lenard- Jones 
potential which realizes the hard-core repulsion, the radius 
of disks R, and the initial density of nanospheres p = 
S/(NttR 2 ) (where S is the simulation area and N is the 
number of considered disks). Varying these parameters in 
reasonable limits all experimentally obtained patterns can 
be successfully modeled [T2] . 

An immediate question that arises in connection with 
the approach based on the spring-block model concerns 
it's computational efficiency. What do we gain by using 
this simplified force pattern instead of the more complete 
variation sketched on Figure [3]? The answer comes from 
several sights. First, the exact shape of the force given 
in Figure [3] is undetermined since the two prefactors for 
the Coulombian and capillarity forces are unknown. Due 
to this, there would appear an extra unknown parameter 
which determines the shape of the resultant force, making 
the model more complex. Secondly, the breaking thresh- 
old imposed on the springs allows the use of a natural 
cutoff for the interparticle forces, and makes simulations 
faster. Third, the simple harmonic force allows a more 
rapid calculation of the resultant forces. Finally, the geo- 
metric condition imposed on springs (not to intersect an- 
other nanosphere) insures the screening effect imposed by 
other nanospheres on capillarity and Coulombian forces. 
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Fig. 4. Characteristic dried nanosphere pattern obtained ex- 
perimentally after the drying front has moved in the horizontal 
direction. 




Fig. 5. The three different regions in the improved model. 

The above mentioned model considered a uniform de- 
crease of the liquid level in the whole system. Experimen- 
tal studies show however that this approach is not valid in 
real systems. During the drying process large liquid level 
gradients are present, leading to the flow of the nanosphere 
suspension from one part of the sample to another one [1, 
[2] . In most of the experimental samples thus drying fronts 
can be observed separating already dried and wet parts 
of the sample. This no n- uniform liquid level is also facili- 
tated by almost all type of suspension deposition methods 
on the silica substrate. The obtained larger-scale patterns 
are strongly influenced by the drying front. Experiments 
have shown that the characteristic larger-scale fracture- 
lines tend to be perpendicular on the direction of the dry- 
ing front (Figured)). 

Recently a simple lattice-gas type model was proposed 
in order to explain the self-assembly process driven by the 
drying liquid front [24J. Monte Carlo simulations with the 
classical Metropolis algorithm proved that patterns resem- 
bling the experimentally obtained ones can be obtained by 
using this elegant approach. 

Our model can be also further improved and made thus 
more realistic to incorporate the effects generated by the 
moving drying front. In order to do this we consider three 
different regions inside the simulated area (Figure [SJ). 

In region I. the liquid level is greater than the na- 
nosphere diameter, no capillary effects are present. This 
models the colloidal suspension state on the substrate (re- 
gions with excess of liquid). In this region the dynamics 




Fig. 6. Snapshot from simulation. On the left side one can 
observe the spheres in region I. and spheres in the right side of 
the picture are already in region III. The drying front is also 
clearly detectable. 

is governed by a hard-core repulsion between spheres and 
a constant force from left to right which is resulting from 
the liquid convection (liquid is flowing from left to right). 
The pattern formation mechanism takes place in region 
II, where the liquid level is lower than the height of the 
nanospheres. This is the region that was considered in our 
earlier model. Here, the dynamics is governed by capillary 
effects, friction with the substrate and repulsion between 
disks, as described earlier. When nanospheres get in this 
region, they become connected with springs having elas- 
tic constants fc^. Finally, the third region (III.) repre- 
sents the final packed monolayer structure where no more 
forces are acting on nanospheres. The separation line be- 
tween region II. and III. defines the drying front in the 
system. 

In the simulations initially we consider all nanospheres 
in region I, and distribute them randomly with a somehow 
smaller density. We define a given thickness for region II 
(usually 1/3 of the width of the simulated region), and 
this region slowly progresses with a constant speed from 
right to the left. Region I gets thus less and less narrow 
and finally we end up with all spheres in region III. Dur- 
ing simulation new particles are constantly supplied from 
the left and driven by the constant force in the direction 
of the separation line between region I and region II. A 
snapshot from the simulation can be observed in Figure 
[6] and a characteristic final pattern is shown in Figure [3 
A qualitative comparison between Figure [4] and Figure [71 
will convince us that this improved model works excel- 
lently and all type of defects observable in experiments 
are reproduced by simulation. 



3 Effect of a random force 

An extra shaking force acting on the nanospheres can be 
now imposed during the simulated drying process. Two 
type of simulations were performed, using both the orig- 
inal and the improved version of the spring-block model. 
In the original version of the spring-block model we con- 
sider an extra stochastic force acting on each nanosphere. 
This random force will model the influence of a random 
shaking applied on the drying nanosphere suspension, and 
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Fig. 7. Typical structure obtained by computer simulations 
with the improved (three region) model. Simulation parameters 
are: F f = 0.003, F b = 0.05, k ini = 0.01, dk = 0.001, the speed 
of the interface between region I and region II was chosen to 
be 10 -6 * L/step (L is the total length of the simulated region 
and steps stands for relaxation steps in the simulation). 

it approximates reasonably well the experimental condi- 
tions of Sasaki and Hanne [9], where an ultrasound bath 
was used to generate the shaking. Second, we consider the 
improved model where in region II. an oscillating force is 
applied on each nanosphere. This force is the same for all 
nanospheres and acts in the direction of the drying front. 
This simulation corresponds to the experimental setup of 
Schope [10]. In this experiment the inclined deposition 
method was used for the deposition of the liquid film, and 
an extra oscillating electric field perpendicular on the di- 
rection of the drag was applied on the nanospheres. 

3.1 The original model with an extra stochastic force 




An extra F r force completely uncorrelated in time and 
space was imposed on each nanosphere. The orientation 
of this force was chosen to be totally random in the sim- 
ulation plane and it's strength was chosen with a uniform 
distribution in a fixed F r G [0, Fr] interval. The value of 
Fr characterizes thus the strength of the applied random 
perturbation. In order to achieve a final stable configura- 
tion it is evident that the condition Fr < Ff must be sat- 
isfied. Fixing the other parameters of the simulations the 
influence of Fr was systematically studied. Fixed bound- 
ary conditions in a disk-like geometry was used and a cen- 
tral rectangular region was studied. Extensive simulations 
proved that by increasing the value of Fr (Fr < Ff ) more 
and more ordered structures are self-assembled. On Figure 
[8] we illustrate this trend for two different parameter sets. 
Simulation results confirm thus that a random shaking 
during the drying process could be helpful in engineering 
more ordered structures. The effect is reasonable, since 



Fig. 8. Influence of the strength of the applied random force 
on the final stable patterns. Results for two different parameter 
sets are shown. The pictures on the left side have been obtained 
with Fb = 0.05 and the ones on the right side with Fb = 0.1. 
Other parameters of both simulations are: Ff = 0.005, kini = 
0.01, dk = 0.001 and p = 0.749. 



the random force will reorder regions during the whole 
simulation process. Dislocations, voids or fracture lines 
formed in the early stages when many springs are present 
will disappear during the dynamics, and the structure will 
evolve more homogeneously. It is also observable from Fig- 
ure [8] that although larger ordered domains are generated 
with increasing Fr values, the fracture lines becomes also 
wider. This effect is however mainly due to the imposed 
fixed boundary conditions. 



6 



F. Jarai-Szabo et al.: Shake induced order in nanosphere systems 




f / F f = 0.5 



F /F f = 0.75 



Fig. 9. Final patterns obtained with an extra harmonic force 
in region II. The successive pictures are for Fo = 0, F//4, F//2 
and 3F//4 amplitude values, respectively. The other simulation 
parameters are as in Figure [8] 



3.2 The improved model with an extra oscillating force 

Within the improved model we applied an additional har- 
monic Fh = Fosin(ut) force in region II. The period of 
the harmonic force was chosen to be of the order of hun- 
dreds of relaxation steps. The Fo amplitude of this force 
was systematically varied in the Fo G [0, Ff) interval. The 
obtained results are similar with the ones obtained with a 
random force. Increasing the value of Fo in the [0, Ff) in- 
terval, the obtained structures become more ordered and 
presents less defects. The result is stable and the quali- 
tative behavior does not change whenever the parameters 
of the simulations are varied. A characteristic series of fi- 
nal patterns as a function of Fo is illustrated on Figure 
The simulated series from Figure [9] describes well the 
experimental results presented in [10J. Simulations sug- 
gested also that higher frequencies eliminate more defects, 
although the qualitative changes in the final structures are 
much less evident than the results presented in Figure [U 
For a systematic study on the frequency dependence sim- 
ulations with much bigger system sizes will be necessary. 



4 Experiments 

As an experimental exercise we have also reproduced the 
results of Sasaki and Hane [9J. Following their work the 
shaking was realized by an ultrasound radiation. The ex- 
perimental samples were prepared following the drop-coat 



method [3,4,5J. Of critical importance for nanosphere or- 
dering is the initial chemical treatment of the glass sub- 
strate in order to render the surface hydrophilic and im- 
prove its wettability (for details see [12J). Polystyrene nanospheres 
of 400 nm diameter, exhibiting negatively charged carboxyl- 
terminated surface with a strongly hydrophobic nature, 
were supplied as monodispersed suspensions in deionized 
water (wt 4%). The original suspension of polystyrene 
nanospheres was diluted by 10 and a volume of 100/iZ di- 
luted solution was evenly spread on the pre-treated sub- 
strates. As the water evaporates and the samples are dried 
the nanospheres self-assemble into close-packed monolayer 
arrays exhibiting many fracture lines, dislocations or other 
type of defects. 

In order to investigate the influence of the shaking ef- 
fect drying was realized in two different conditions. Six 
samples were studied, half of them were dried simply in 
the normal atmosphere of the lab while the other half was 
exposed to a shaking effect induced by an ultrasound bath. 

We have used an Elma Transsonic 35 kHz frequency 
ultrasound bath filled with water. The glass plates holding 
the nanosphere solution was positioned on the surface of 
the water by using a 10 cm diameter and 0.4 cm thick 
cork floater with a disk-like geometry. In the middle of the 
floater a hole with a 0.8 cm diameter was created, allowing 
direct contact between the liquid from the ultrasonic bath 
and the surface of the glass plates. Shaking is induced thus 
by the acoustic vibrations transmitted through the glass 
plates. 

In both cases (with and without ultrasound) the sam- 
ples were completely dried. Drying was achieved in ap- 
proximately 45 minutes for the sonicated samples. For 
the nonsonicated samples complete drying was achieved 
roughly two times slower. The microstructure of the dried 
samples were studied by scanning electron microscopy (SEM) 
using a JEOL JSM 5600 LV electronic microscope. By sur- 
veying different regions of the substrates we have qualita- 
tively analyzed and compared the engineered nanosphere 
patterns. 

Our experimental exercise confirms the picture sug- 
gested by Sasaki and Hane: sonicated samples shows more 
extended and ordered triangular lattice structures than 
non-sonicated ones. From a first look on a randomly se- 
lected region on a sonicated and a non-sonicated sample, 
this result is however not immediate. One has to remember 
that due to the non-homogeneous spreading of the liquid 
film, different regions in the same sample might exhibit 
different level of ordering. In order to obtain thus a clear 
conclusion many different regions have to be analyzed and 
a statistical conclusion has to be drawn. We have analyzed 
tens of different regions and qualitatively found that soni- 
cated samples had less number of fracture and dislocation 
line, the triangular lattice structures being more extended 
in this case. For visually illustrating this, on Figure fT()l we 
present a characteristic pattern found on sonicated sam- 
ples in comparison with a characteristic pattern found on 
non- sonicated samples. 
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Fig. 10. Visual comparison between characteristic regions of 
non-sonicated (left) and sonicated samples (right). 



As an exercise we have also repeated the experiments 
done by Sasaki and Hane, confirming their results. Dur- 
ing the experiments we realized that in order to get more 
extended triangular lattice structures and less defects one 
will have to further optimize the experimental conditions. 
To make the experiments less costly computer simulations 
could be useful. The present work shows that a simple 
Burridge-Knopoff type model works empirically well and 
can be helpful to give a first estimate on the influence of 
many experimentally controllable parameters. 
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5 Conclusions 

Nanospehere self-assembly from a drying suspension was 
studied in the presence of a moderate intensty shaking. To 
understand this pattern formation phenomenon two mod- 
els were considered, (i) the simple Burridge-Knopoff type 
spring-block model [12] and (ii) an improved version of 
this model where the withdrawal of the liquid front and 
it's effect on self-assembly is also taken into account. In 
both cases beside the relevant capillarity, Coulombian and 
friction (pinning) forces and also an extra random shak- 
ing was imposed. Computer simulations suggested a clear 
conclusion: a moderate intensity extra shaking imposed on 
the drying nanosphere system will yield more ordered final 
patterns. Increasing the intensity of the random force up 
to the limit of the pinning forces acting on the nanospheres 
is always beneficial and eliminates most of the defects. The 
simulation performed on our original spring-block model 
with an extra random force describes well the experimen- 
tal results of Sasaki and Hane [9], while the simulations 
realized on the improved model with an extra directional 
harmonic force acting in region II corresponds to the ex- 
perimental setup of Schope [10]. In agreement with ex- 
periments computer simulations also yield that primary 
(larger) fracture lines are formed predominantly in the 
direction of the moving liquid front (for experiments see 
Figure [4] and for simulation Figure [6] or [7)) . The patterns 
and the trend obtained by computer simulations describes 
thus qualitatively well the experimentally obtained ones, 
giving new evidences for the applicability of the simple 
Burridge-Knopoff type models in modeling this system. 

We have not analyzed here the influence of the fre- 
quency of the shaking. The first reason for this is that 
the influence of the shaking frequency on the final pat- 
tern is less evident than the influence of the intensity of 
the shaking force. In order to get reliable data computer 
simulations on much larger systems have to be done. The 
computing power presently available for us does not al- 
low studying much larger systems in reasonable (months) 
computing time. Second, the simulations in the present 
form lack real time. Instead of time we used the relax- 
ation steps for characterizing the dynamics. Relaxation 
steps might correspond however to different time intervals 
and in this sense the frequency of the imposed oscillating 
force has only a qualitative meaning. 



The present work has been supported by the research 
grant CNCSIS 41/183 and CEEX-Nanobiospec 2006. The 
research of F. Jarai-Szabo has also been supported by the 
Felowship Program for Transborder Hungarian Scientific 
Research - Hungarian Academy of Science. Discussions 
and helpful comments from Prof. Peter Kralchevsky are 
gratefully acknowledged. 



References 

1. N.D. Denkov, O.D. Velev, P.A. Kralchevsky, LB. Ivanov, H. 
Yoshimura and K. Nagayama, Langmuir 8, 3183 (1992) 

2. N.D. Denkov, O.D. Velev, P.A. Kralchevsky, LB. Ivanov, H. 
Yoshimura and K. Nagayama, Nature 361, 26 (1993) 

3. Ch. L. Haynes, R.P. van Duyne, J. Phys. Chem. B 105, 
5599 (2001) 

4. K. Kempa et al., NanoLetters 3, 13 (2003) 

5. W.A. Murray, S. Astilean, W.L. Barnes, Phys. Rev. B 69, 
165407 (2004) 

6. A. A. Chabanov, Y. Jun. D.J. Norris, Appl. Phys. Lett. 84, 
3573 (2004) 

7. E. Vasco Appl. Phys Lett. 85, 3714 (2004) 

8. K. Chen, A. Taflove, Y. L. Kim and V. Backman, Appl. 
Phys. Lett. 86, 033101 (2005) 

9. M. Sasaki and K. Hane, J. Appl. Phys. 80, 5427 (1996) 

10. H.J. Schope, J. Phys., Condens. Matter 15, L533 (2003) 

11. M. Baia, L. Baia and S. Astilean, Appl. Phys. Lett. 88, 
143121 (2006) 

12. F. Jarai-Szabo, S. Astilean and Z. Neda, Chemical Physics 
Letters 408, 241 (2005) 

13. R. Burridge and L. Knopoff, Bull. Seis. Soc. Amer. 57, 341 
(1967) 

14. J. H. E. Cartwright, E. Hernandez-Garcia and O. Piro, 
Phys. Rev. Lett. 79, 527 (1997) 

15. K.-t. Leung and Z. Neda, Phys. Rev. Lett. 85, 662 (2000) 

16. Z. Neda, K.-t. Leung, L. Jozsa and M. Ravasz, Phys. Rev. 
Lett. 88, 095502 (2002) 

17. D.Y. Chan, J.D. Henry and L.R. White, J. Colloid Inter- 
face Sci. 79, 410 (1981) 

18. P.A. Kralchevsky, V.N. Paunov, N.D. Denkov, LB. Ivanov 
and K. Nagayama, J. Colloid Interface Sci. 155, 420 (1993) 

19. V.N. Paunov, P.A. Kralchevsky, N.D. Denkov and K. Na- 
gayama, K. Nagayama, J. Colloid Interface Sci. 157, 100 
(1993) 

20. CD. Dushkin, P.A. Kralchevsky, H. Yoshimura and K. Na- 
gayama, Phys. Rev. Lett. 75, 3454 (1995) 



8 



F. Jarai-Szabo et al.: Shake induced order in nanosphere systems 



21. CD. Dushkin, P.A. Kralchevsky, V.N. Paunov, H. 
Yoshimura and K. Nagayama, Langmuir 12, 641 (1996) 

22. CD. Dushkin, H. Yoshimura and K. Nagayama, J. Colloid 
Interface Sci. 181, 657 (1996) 

23. K.D. Danov, B. Pouligny, and P.A. Kralchevsky, Langmuir 
17, 6599 (2001) 

24. E. Rabani, D.R. Reichman, P.L. Geissler and L.E.Brus, 
Nature 426, 271 (2003) 



